home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / clacon.z / clacon
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. CCCCLLLLAAAACCCCOOOONNNN((((3333FFFF))))                                                          CCCCLLLLAAAACCCCOOOONNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CLACON - estimate the 1-norm of a square, complex matrix A
  10.  
  11. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  12.      SUBROUTINE CLACON( N, V, X, EST, KASE )
  13.  
  14.          INTEGER        KASE, N
  15.  
  16.          REAL           EST
  17.  
  18.          COMPLEX        V( N ), X( N )
  19.  
  20. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  21.      CLACON estimates the 1-norm of a square, complex matrix A.  Reverse
  22.      communication is used for evaluating matrix-vector products.
  23.  
  24.  
  25. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  26.      N      (input) INTEGER
  27.             The order of the matrix.  N >= 1.
  28.  
  29.      V      (workspace) COMPLEX array, dimension (N)
  30.             On the final return, V = A*W,  where  EST = norm(V)/norm(W) (W is
  31.             not returned).
  32.  
  33.      X      (input/output) COMPLEX array, dimension (N)
  34.             On an intermediate return, X should be overwritten by A * X,   if
  35.             KASE=1, A' * X,  if KASE=2, where A' is the conjugate transpose of
  36.             A, and CLACON must be re-called with all the other parameters
  37.             unchanged.
  38.  
  39.      EST    (output) REAL
  40.             An estimate (a lower bound) for norm(A).
  41.  
  42.      KASE   (input/output) INTEGER
  43.             On the initial call to CLACON, KASE should be 0.  On an
  44.             intermediate return, KASE will be 1 or 2, indicating whether X
  45.             should be overwritten by A * X  or A' * X.  On the final return
  46.             from CLACON, KASE will again be 0.
  47.  
  48. FFFFUUUURRRRTTTTHHHHEEEERRRR DDDDEEEETTTTAAAAIIIILLLLSSSS
  49.      Contributed by Nick Higham, University of Manchester.
  50.      Originally named CONEST, dated March 16, 1988.
  51.  
  52.      Reference: N.J. Higham, "FORTRAN codes for estimating the one-norm of a
  53.      real or complex matrix, with applications to condition estimation", ACM
  54.      Trans. Math. Soft., vol. 14, no. 4, pp. 381-396, December 1988.
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.